VARIATIONAL FORMULATION FOR THE KPZ AND RELATED KINETIC 

EQUATIONS 
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We present a variational formulation for the Kardar-Parisi-Zhang (KPZ) equation that leads to a 
thermodynamic-like potential for the KPZ as well as for other related kinetic equations. For the 
KPZ case, with the knowledge of such a potential we prove some global shift invariance proper- 
ties previously conjectured by other authors. We also show a few results about the form of the 
stationary probability distribution function for arbitrary dimensions. The procedure used for KPZ 
was extended in order to derive more general forms of such a functional leading to other nonlinear 
kinetic equations, as well as cases with density dependent surface tension. 
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1. Introduction 

Phenomena far from equilibrium are ubiquitous in na- 
ture, including among many other, turbulence in fluids, 
interface and growth problems, chemical reactions, bi- 
ological systems, as well as economical and sociological 
structures. During the last decades the focus on sta- 
tistical physics research has shifted towards the study 
of such systems. Among those studies, the understand- 
ing of growing kinetics at a microscopic as well as on a 
mesoscopic level constitutes a major challenge in physics 
and material science [Tong & Williams, 1994; Barabasi 
& Stanley, 1995; Halpin-Healy & Zhang, 1995; Marsili 
et ai, 1996]. Some recent papers have shown how the 
methods and know-how from static critical phenomena 
have been exploited within nonequilibrium phenomena of 
growing interfaces, obtaining scaling properties, symme- 
tries, morphology of pattern formation in a driven state, 
etc [Hentschel, 1994; Prahofer & Spohn, 2004; Lopez et 
ai, 2005; Fogedby, 2006; Ma et ai, 2007; Castro et ai, 
2007]. 

Even though it was (briefly) discussed in [Cross & Ho- 
henberg, 1993], there is still a common belief that the 
nontrivial spatial-temporal behavior occurring in sev- 
eral nonequilibrium systems, originates from the non- 
potential (or non-variational) character of the dynam- 
ics, meaning that there is no Lyapunov functional for 
the dynamics. However, Graham and co-workers have 
shown in a series of papers [Graham, 1987; Graham & 
Tel, 1990] that a Lyapunov-like functional exists for a 
very general dynamical system, the complex Ginzburg- 
Landau equation. Such a functional is formally defined 
as the solution of a Hamilton-Jacobi-like equation, or ob- 
tained in a small gradient expansion [Descalzi & Graham, 
1994]. The confusion associated with the qualification of 
nonvariational dynamics comes from the idea that the 
dynamics of systems having nontrivial attractors (limit 
cycle, chaotic) cannot be deduced from the minimization 
of a potential playing the same role as the free energy 
in equilibrium systems [Cross & Hohenberg, 1993]. Nev- 
ertheless, this does not preclude the existence of a Lya- 
punov functional for the dynamics that will have local 
minima identifying the attractors of the system [Mon- 



tagne et ai, 1996; Wio, 1997; Wio et ai, 2002; Wio & 
Deza, 2007]. However, once the system has reached an 
attractor that is not a fixed point, the dynamics pro- 
ceeds inside the attractor driven by nonvariational con- 
tributions to the dynamical flow, that do not change the 
value of the Lyapunov functional, implying that the dy- 
namics is not completely determined once the indicated 
functional is known. This situation has known exam- 
ples even in equilibrium statistical mechanics [Hohen- 
berg & Halperin, 1977]. Hence, the Lyapunov functional, 
or nonequilibrium potential (NEP) [Graham, 1987; Gra- 
ham & Tel, 1990; Wio, 1997; Ao, 2004], plays the role 
in nonequilibrium situations of a thermodynamical-likc 
potential characterizing the global properties of the dy- 
namics: attractors, relative (or nonlinear) stability of 
these attractors, height of the barriers separating attrac- 
tions basins, offering the possibility of studying transi- 
tions among the attractors due to the effect of (thermal) 
fluctuations. 

In a recent series of papers we have shown several results 
related to the obtention of the indicated NEP's for other 
system's classes: scalar and non-scalar reaction-diffusion 
systems [Bouzat & Wio, 1999; von Haeften et at, 2004; 
von Haeften et at, 2005; von Haeften & Wio, 2007]. In 
particular we have exploited those results for the study 
of stochastic resonance [Gammaitoni et ai, 1998] in ex- 
tended systems [Wio, 1997; Bouzat & Wio, 1999; von 
Haeften et ai, 2000; von Haeften et ai, 2004; von Haeften 
et ai, 2005; von Haeften & Wio, 2007; Tessone & Wio, 
2007; Wio & Deza, 2007;]. In those works, we have an- 
alyzed problems of stochastic resonance in scalar and 
activator-inhibitor systems, in systems with local and 
nonlocal interactions, studied system-size stochastic res- 
onance, etc. 

Here, and related to the kinetics of growing interfaces, 
we discuss the case of the Kardar-Parisi-Zhang equa- 
tion (KPZ) [Kardar, Parisi & Zhang, 1986; Medina et 
ai, 1989]. This equation, that describes the evolution of 
h(x, t), a field that corresponds to the height of a fluctu- 
ating interface, reads 



dh(x,t) 
dt 



vV 2 h{x,t) 



2 



where £(x, t) is a Gaussian white noise, of zero mean 
£)) = 0) and correlation (£(x, t)f (as', t')) = 2e^(x - 
x')6(t — t'). As indicated above, this nonlinear differen- 
tial equation describes fluctuations of a growing interface 
with a surface tension given by is, A is proportional to the 
average growth velocity and arises because the surface 
slope is parallel transported in such a growth process. 
Opposing to a claim in a recent paper [Fogedby, 2006]: 
The KPZ equation is in fact a genuine kinetic equation 
describing a nonequilibrium process in the sense that the 
drift v\7 2 h + -|V/i • V/i — F cannot be derived from an 
effective free energy; in [Wio, 2007] it was shown that 
such a nonequilibrium thermodynamic-like poten- 
tial (NETLP) for the KPZ equation exists. 
In this paper we present the approach to obtain the 
NETLP for the KPZ equation, and also show how it is 
possible to obtain such a NETLP for other related ki- 
netic equations, as well as extend the procedure to other 
general situations. The organization is as follows. In the 
next Section we show how to obtain the NETLP for the 
KPZ case and, exploiting its knowledge, we discuss con- 
jectures advanced in [Hentschel, 1994] and how they are 
fulfilled. Next we discuss about the form of the stationary 
(or asymptotic) probability distribution function, valid 
for any dimension. In the following Section we discuss 
the situation with a non-local interaction and symmetric 
kernels, showing how this form of extending the deriva- 
tion procedure allows to obtain more general forms of 
kinetic equations. The following Section discusses the 
case of density dependent surface tension, while the last 
one presents some final comments. 

2. Kardar-Parisi-Zhang case 
2.1 Derivation 

In order to show how to obtain the above indicated func- 
tional, we start considering the following simple (gen- 
eral) scalar reaction-diffusion equation with multiplica- 
tive noise 

J^(x, t) = vV 2 <}>(x, t) + f(cj>{x, t)) + <f>(x, t)n(x, t), (2) 

where f(cj>(x,t)) is a general nonlinear function. The 
noise n(x,t) is a Gaussian white noise, of zero mean 
and intensity a. For this equation we assume the 
Stratonovich interpretation. 

It is known that the system in Eq. ([2]) has the following 
NETLP 



(x,t) 



/foOdp+£(V#x,t)) a \dx, (3) 



in I Jo 2 

where indicates the integration range, and 

J^(M) = -^^ + #Mh(M); (4) 



usual. As has been shown in previous works [Izus et al, 
1998; Bouzat & Wio, 1998], it also fulfills the Lyapunov 
characteristic < (it is worth here commenting 

that, as a matter of fact, this conditions is only valid in 
a weak noise limit). 

Exploiting the so called Hopf-Cole transformation we 
now define a new field, h(x,t) that, as indicated before, 
corresponds to an interface height, 



h(x, t) 



2v 



In 4>(x, t), 



(5) 



with the inverse 



(j>(x,t) 



As 4>(x,t) > 0, h(x,t) is always well defined. The trans- 
formed equation reads 



A 



A 



d t h(x,t) = v^h+- (Vhy + —e--f(h) + ax,t), (6) 

where f(h) = f(<p) = f(e^ h ). Now, in order to reduce 
our general result to the one shown in Eq. ([6]) , we assume 
f(<f>(x, t)) = a<j)(x, t) (with a some constant). In this way 
Eq. ([2]) becomes Eq. (fT]), with a = jr;K and a = 

(^j) 2 £■ However, the noise term that in the original 
equation (Eq. @) has a multiplicative character, in the 
transformed equation (Eq. (fT])) becomes additive. 
If we now apply the same transformation to the NETLP 
indicated in Eq.Q, restricting ourselves to f(<p(x,t)) — 
a<p(x,t), we obtain 



T[h} = 



)h{x,t) 



£2 
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~(Vh(x,t)) 2 



dx. (7) 



It is easy to prove that this functional fulfills both, the 
relation 



5^" " - r ^MB) +m ' } '' 



(S) 



as well as the Lyapunov characteristic -^!F[h\ < 0, where 
the function T[h] is given by 



T[h] 



Hence we have a free energy-like functional from 
where the KPZ kinetic equation can be obtained through 
functional derivation. Clearly, the contribution to the 
variation coming from the boundaries is again null. 
It is worth to consider once more Eq. ^ , but now assum- 
ing f(4>(x, t)) = a<j)(x, t) — b(j){x, t) 3 , i.e. we include a lim- 
iting (or saturating) term. In such a case, the reaction- 
diffusion equation corresponds to a form of the so called 
Schlogl model [Mikhailov, 1990; Wio, 1994]. The "po- 
tential" contribution to the NETLP in Eq. © has the 
form 



where the contribution from the boundaries is null, due to 
the variation 5<p being fixed (= 0) at these boundaries, as 
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Applying again the indicated Hopf-Cole transformation, 
in Eq. ([T]) a new associated term arises, having the form 
— jei h ( x,t \ with b — ^7. The new kinetic equation cor- 
responds to a form of the so called bounded-KPZ [Grin- 
stein et ah, 1996; Tu et at, 1997; de los Santos et at, 
2007]. Clearly, we will also have an extra term in the 
associated NETLP (Eq. ([7])). However, here we restrict 
ourselves to the case 6 = 0, and only analyze the more 
"usual" form of the KPZ equation indicated by Eq. ([TJ) . 

2.2 Some Properties 

Let us now check some of the properties assumed previ- 
ously for such a functional. According to the analysis of 
global shift invariance in [Hentschel, 1994], it is easy to 
see that the relations indicated by Eq. (9) in the indi- 
cated paper are fulfilled. That is, we can prove that if I 
is an arbitrary (constant) shift 



F[h + l] = K[l]T[h] 
T[h + l] = Kllj-^lh], 



(9) 



with K[l] = e^ 1 = (^-YT[l]-\ 

To prove other conjectures also indicated in [Hentschel, 
1994], we introduce the free energy-like density !F[h, V/i], 
defined through 



T[h] = J dxT[h,Vh]. 



Hence, T[h,Vh] = ±ei h( - s ^ [-K + ^(Vh(x, t)) 2 ]. The 
relations we refer are 



f[h, V/i] 

r[h,vh] 



sh 



Ti[(V/i) 2 



-sh 



(10) 



and according to the form of J-[h, V/i], it is clear that 
the first relation above results obviously true, while for 
the second relation we have that T[h, V/i] = e - *'*( a! > t )r o , 
where r o = 1, and s — as T[h] is not function of V/i. In 
addition, it can be also proved that the indicated NETLP 
is also invariant under the nonlinear Galilei transforma- 
tion that, as discussed in [Fogedby, 2006], are fulfilled by 
the KPZ equation. 

It is here adequate to make a warning. We have found 
that, from the indicated free energy-like functional for the 
KPZ kinetic equation and by a functional derivative, we 
can obtain a form that resembles a (relaxation) model A 
according to the classification in [Hohenberg & Halperin, 
1977]. In one hand, in the standard "model A" it is 
known that the dynamics can be seen as a superposition 
of modes that decay exponentially towards a steady state. 
In addition, it is also known that the time-dependent cor- 
relations obey certain constraints such as positivity. On 
the other hand, in the KPZ problem it is known that 
the relaxation of perturbations decay in a stretched ex- 
ponential manner [Schwartz & Edwards, 2002; Edwards 
& Schwartz, 2002; Colaiori & Moore, 2002; Katzav & 



Schwartz, 2004]. Hence, even though Eq. (|8|) looks simi- 
lar to "model A" , its behavior is far from trivial, and we 
have no a-priori intuition to what its dynamic could be. 
Clearly, this is a point to have in mind when suggesting 
any kind of Antsatz for the temporal behavior. 

3. Probability Distribution Function for KPZ 
3.1 General Aspects 

The knowledge of the NETLP indicated in the previous 
section, allows us to readily write the asymptotic long 
time probability distribution (pdf). We start writing the 
(functional) Fokker-Planck equation (FPE) associated to 
Eq. (TT]), that reads 

^P[h(x,t)] =^eteJ-{-^V 2 /i(x,i) 

+ ^ (Vh(x, t)) 2 ]v[h(x, *)] + e ^V[h(x, *)]}; (11) 

where, in order to focus in the most relevant aspects and 
to simplify, we adopted K = 0. 

In Chap. 6 of [Barabasi & Stanley, 1995] as well as in 
Sect. 3.5 of [Halpin-Healy & Zhang, 1995] the form of 
the probability distribution function (pdf) for the one di- 
mensional KPZ equation was discussed. The form of such 
a pdf for an arbitrary dimension, so far, is not known. 
However, here we show the general form of such a pdf. 
The knowledge of such a NETLP for the KPZ equation 
allows us to readily write the asymptotic long time prob- 
ability distribution function (pdf), valid for (worth to be 
remarked) any dimension, which (due to the "diago- 
nal" character of T[h]) is given by 



V as [h{x)} 



-{-;/*£> rM ^} 

cxp|-iy dx (v[h]T[h] 



m , , ST[ip] 



T 



cxpj — ^- J dx(Vh) 2 



A f , f Hs) 
-— ax 
2eJ J href 



d^(V^) 2 } 



exp 



(12) 



where, we reiterate, we assumed K a = 0, and h re f is an 
arbitrary reference state. The second line results by us- 
ing functional methods (see for instance [Hanggi, 1985]). 
The third line shows a nice structure, where we can iden- 
tify a contribution, with a Gaussian dependence on the 
slope, plus a "correction" term proportional to A. The 
above indicated result is the valid (albeit "formal") so- 
lution for arbitrary dimension irrespective of boundary 
conditions, while for the one-dimensional case, and the 
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adequate (periodic) boundary conditions, it is possible to 
show that this pdf reduces to the well known Gaussian 
result [Barabasi & Stanley, 1995; Halpin-Healy & Zhang, 
1995]. 

Clearly, a relevant point is related to the interpretation of 
the integral over the function ip(x) in Eq. (fT2|) . In order 
to simplify we consider the 1-d case, and focus only on 
the first term. Using functional techniques, we have 




2. \dx 



where we have used that 5 ^p^J ) — S(x — x'). Another 
way to grasp it is via a discrete representation 



-Y 

2 ^ 



#7 TT Y t^+i ~ ^] 



that, as only the / = j and I = j — 1 terms survive, yields 

« ? E fo+i - ( 14 ) 



In both cases we obtain the known result. 
The interpretation for the second (KPZ) term is analo- 
gous to the one in the example indicated above. However, 
for this case the situation is a little more delicate. Several 
recent papers have discussed different alternatives to the 
discrete form of such a contribution in the KPZ equa- 
tion in 1-d [Newman & Bray, 1996; Lam & Shin, 1998; 
Ma, Jiang & Yang, 2007]. We will not come into these 
details here, that will be deeply discussed in [Revelli & 
Wio, 2008]. 

3.2 Nonequilibrium Potential for KPZ 

The last line in Eq. (fT!?]) defines the functional 



V 



dx (V/i) 2 - - dx dip (V^) 2 . (15) 
2 ./ 2 J J Ksf 



The variation of this functional give us 

d , , s 6$[h] „ , 



(16) 



This functional also fulfills the Lyapunov condition 

Hence, such a functional could be identified as the 
nonequilibrium potential [Graham, 1987; Graham & Tel, 



1990; Wio et ai, 2002] for the KPZ case. More, even 
though it could be also identified with a Hamiltonian for 
KPZ, it is worth to remark that it differs from the form 
indicated in Eq. (3.4) of [Halpin-Healy & Zhang, 1995]. 
The last functional form (albeit so far only formal) could, 
in principle, allow us to exploit several known techniques 
[Graham, 1987; Graham & Tel, 1990; Wio et at, 2002]. It 
also shows that the claim by previous authors indicated 
at the Introduction is not true. 

3.3 NEP's application: A simple example 

Here we discuss a simple example in order to show the 
possibilities that offers the knowledge of such a NETLP. 
For this example we analyze a sightly different situation 
than the one studied in Eq. ([T]) and its associated NEP 
in Eq. (|15p . We only consider a spatial quenched noise 
( "disorder" ) instead of the spatial-temporal noise consid- 
ered so far. The equation associated to such a problem 
is 

^-h{x,t) = vV 2 h{x,t) + ^{Vh{x,t)f 



2 



(17) 



where, as in previous studies [Nattermann & Renz, 1989; 
Krug & Halpin-Healy, 1993; Ramasco, Lopez & Ro- 
driguez, 2006; Szendro, Lopez & Rodriguez, 2007], d(x) 
is a quenched, Gaussian distributed, noise. For this case 
we have that 



d_ 

di 



h{x,t) = -T[h] 



5T[h] 
Sh(x, t) 



(18) 



with the same form of J-\h] as in Eq. (JT]), but K replaced 
byK + §{x). 

The indicated studies have shown that the front profile 
presents a triangular structure [Ramasco, Lopez & Ro- 
driguez, 2006; Szendro, Lopez & Rodriguez, 2007] and, 
clearly, it is of relevance to determine the slope of such 
structures. Using the known form of the NETLP, we can 
minimize this free-energy-like functional and obtain, in 

1 /3 

the 1-d case, that such a slope is a = (jQ) , a value 
that agrees quite well with the numerical evaluations. 

4. Other Kinetic Equations: Non Locality 

We can go still further and look for the possibility of de- 
riving a NETLP for more general forms of kinetic equa- 
tions. To this end, let us assume that we have the fol- 
lowing non-local reaction-diffusion equation [Wio et al., 
2002] 

d 

— <f)(x,t)= v \7 2 (/)(x, t) + a<l)(x, t) 



[3 / dx'G{x,x')(j)(x',t) 
Jn 

+<j)(x,t)r)(x',t), 



(19) 



where, as discussed in [Bouzat & Wio, 1999; von Haeften 
et al, 2004; von Haeften & Wio, 2007], the kernel 



5 



G(x,x') could be of a very general character, and j3 is 
the interaction intensity. It was shown that the form of 
the associated NETLP is 

m= Jj-~<i>(x,t) 2 +~(vct>(x,t)) 2 

+0 / dx' cj>{x,t)G(x,x')4>{x',t)\dx. (20) 

As we have done before, using the Hopf-Cole transforma- 
tion we obtain a generalized-non-local form of the KPZ 
equation 

^h(x, t) = vV 2 h{x, t) + ^ (Vh{x, t)) 2 + K 



-p e -*sKs,*) I dx'G(x,x')e* 
Jn 



(21) 



Even though the nonlocal contribution indicated above 
differs from those discussed in [Mukherji & Bhattachar- 
jee, 1997; Katzav, 2003], it is clear that such form of a 
nonlocal term is also of great interest. Repeating the 
previous procedure we find the associated NETLP 



F[h] 



dxe^^ 



A 

~2v' 



-Kn 



^) (Vfc(M)) a 



+ p e -£ h (z,t) I dx'G(x,x')e* 
In 



h(x',t) 



(22) 



At this stage it is required to make some assumptions 
about the kernel. 

We assume that the nonlocal kernel has translational in- 
variance, that is G(x,x') = G(x — x'). Also, that it is of 
(very) "short" range, that allows us to expand it as 



G(x -x') = ^T A 2n 8 {2n \x - x'), 



(23) 



with 5( n ) (x— x') = V"/5(x— x'), and symmetry properties 
taken into account. Exploiting this form of the kernel, 
we arrive to the following contributions in Eq. (|21[) 



e 2, 



■h(x,t) 



G{x - x')e& h{i '' t} 




-V l h 



2v 



(»[ (Vh) 2 W 2 h 

+ 2(A) 2 v2(v , )2 _(M 2 (v2 , r 



2v 



A. 



(24) 



where the last term indicates contributions of order n > 3 
(2n = 6). The parameter (3 could (in principle) be pos- 
itive or negative, indicating an inhibitor or an activator 
role for the nonlocal interaction term, respectively. 



These contributions have the same form of those ones 
that arose in several previous works, where scaling prop- 
erties, symmetry arguments, etc, have been used to 
discuss the possible contributions to a general form of 
the kinetic equation [Hentschel, 1994; Linz et ai, 2000; 
Lopez et ai, 2005]. Clearly, the different contributions 
that arose in Eq. (f2~4")l are tightly related to several of 
other previously studied equations, like the Kuramoto- 
Sivashinsky [Sivashinsky, 1977; Kuramoto, 1978], the 
Sun-Guo-Grant equation [Sun et at, 1989], and others 
[Hentschel, 1994]. 

5. Density Dependent Surface Tension 

We can also go further in another direction. Let us con- 
sider the case of density dependent diffusion ( surface ten- 
sion). Following the same original approach, we start 
considering the following reaction-diffusion equation with 
a density dependent diffusion for a field <f>(x, t), as studied 
in [von Haeften et ai, 2000; Tessone & Wio, 2007] 

dt^x, t) = V (K4W) + f{4>) + cf>(x, t) V (x, t). (25) 

As before, /((/>) is a general nonlinear function, that 
we collapse to f{4>) — a<f>(x,t). As in Eq. @, r](x,t) 
is a delta correlated white noise of intensity a. We 
assume vicj)) — v Q g{(j)), that -using the Hopf-Cole 
transformation- yields v(h) — v 9{h)- 
The transformed equation reads 

d t h(x,t) = V(i>(/i)V/i) + ^ (Vh) 2 +K +£(x,t), (26) 



with X(h) = \ g(h) and a = -^-K . This equation has 
a KPZ-like form, with both a density dependent surface 
tension term and a density dependent nonlinear coupling, 
with both nonlinearities having the same functional de- 
pendence. 

In order to see the associated NETLP, let us remember 
its form for the original reaction-diffusion equation (Eq. 
JU}), it reads 



m -/„{-jT 



fWW)<#' + ^(K$V$ a | dx, 



(27) 



that fulfills [von Haeften et ai, 2000] 

Applying once more the Hopf-Cole transformation, we 
obtain 



r\i,\ : / dx ( { - / <u,( —i-[it\K„ 



+ \ (i>(/ i )e^V/x) 2 }(28) 
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From this NETLP we can prove that 

= ~ t[h] wB) + ^ (29) 

where £(x,t) is a white noise of intensity e (with a = 
(^7") e )' anc ^ tnc mnc tion f [h] is given by 

'2v n \ e 



T[h] = 



2 -^fe(ii) 



A Q J i>{h) 



As before, if we have f((j>(x, t)) = a<f>(x, t) — b <f>(x, t) 3 , we 
get a generalized form of the "bounded KPZ" equation. 



As indicated in the literature, dynamic renormalization 
group techniques, being useful and powerful, in many 
cases only offers incomplete results, having no access 
to the strong coupling phase [Barabasi & Stanley, 1995; 
Wiese, 1998]. Hence, it is clear the need of alternative 
ways to analyze the KPZ and related problems, as for in- 
stance the self-consistent expansion [Katzav & Schwartz, 
1999; Katzav, 2003]. The present results open new possi- 
bilities of making non-perturbational studies for the KPZ 
problem. For instance, through the analysis of long time 
mean values of h(x,t). In a similar way, it would be 
possible to obtain correlations, and from them to extract 
information about scaling exponents. Such study will be 
the subject of forthcoming work. 



6. Conclusions 

The previous results indicates that, for a very general 
form of a KPZ-like equation, we can device a NETLP a 
la carte. Consider the following equation 

d t h(x, t) = vV 2 h + ^ (V/i) 2 + F(h) + £(x, t) (30) 

where F(h) is a general nonlinear function of h. Ac- 
cording to the previous results we could readily write the 
associate NETLP that has the form 



?[h] = 



A 

2^ 



e " F(u)du 



+ 



81/ 



/ Ah 



V/i 



dx. 



(31) 



Clearly, it fulfills the Lyapunov condition j^T[h{x, t)] < 
0, and also, through functional derivation, allows us to 
obtain the kinetic equation as 



d_ 

dt 



h(x,t) = -T[h] 



6F[h] 
5h(x, t) 



(32) 



with T[h] = (2f) 2 e-^'*). Hence, it is clear that we 
can write the NETLP for a very general KPZ-like form, 
independently of the actual form of F(h). 
Summarizing, we have here found the form of the Lya- 
punov functional or NETLP for the KPZ equation. More, 
we have devised a way to extent the procedure to de- 
rive it, and in such a way we were able to derive more 
general forms, including several kinetic equations studied 
in the literature of interface growing phenomena. Even 
the case of density dependent surface tension, have been 
discussed. From this NETLP, and through a functional 
derivative, we have obtained either the KPZ as well as 
other generalized kinetic equations. We have also shown 
that the NETLP for KPZ fulfills global shift properties, 
as well as other ones anticipated for such an unknown 
functional. More, we have found the exact expression, 
valid for any dimension, for the probability distribution 
function, and have commented on the result of a simple 
example that indicate the usefulness of such a functional. 
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